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ABSTRACT 
The general properties of the gamma distribution, which has several applications in meteorology, are discussed. 


A short review of the general properties of good statistical estimators is given. 
tion to show that the maximum likelihood estimators are jointly sufficient. 


This is applied to the gamma distribu- 
A new, simple approximation of the 


likelihood solutions is given, and the efficiency of the fitting procedure is computed. 


1. INTRODUCTION 


In 1947, the writer [1] developed approximate solutions 
of the maximum likelihood (M.1L.) equations for the incom- 
plete gamma distribution commonly called the gamma 
distribution. The purpose of this note is to give the 
development of the estimates which have since been widely 
employed. The first application of the methods was to 
rainfall data [2]; later it was found that the gamma distri- 
bution has wide application in meteorology to problems 
where the climatological variable has a physical lower 
bound of zero but no nonstatistical upper bound. 

The gamma distribution is a special case of the Pearson 
Type III distribution where the locus parameter is zero. 
Fisher [3] first gave the M. L. equations for this distribu- 
tion; however, as is often the case with M. L. estimation 
the equations are not conveniently solved. Our approxi- 
mations make the estimation of the distribution param- 
eters hardly more difficult than the estimation of the mean 
and standard deviation of the normal distribution. 


2. THE GAMMA DISTRIBUTION AND PROPERTIES 


The gamma distribution is a 2-parameter frequency 
distribution given by the equation 


-2/B. 
f(z) z > (1) 
Here z is the random variable, 8 scales z and is therefore 
the scale parameter, y is the shape parameter, [ is the 
usual gamma function, and f(z)=0 forz<0. It will be 
noted that the distribution has a zero lower bound and is 
unlimited on the right. It is positively skewed, the 
amount of skew depending inversely on the shape factor y. 
The mode of the distribution is at 8 (y—1) if y>1 and 
at zero if 0<y<1. In the latter case, the distribution 
is J-shaped. For y=1 the distribution is exponential 
with ordinate 1/8 at z=0; for y<.1 the ordinate at z=0 
is infinite. The gamma distribution is closely related to 
the chi square distribution, for x?/2 is a gamma variate 
with y= and B=1. 
The moments about zero of the gamma distribution are 
given by the relation 


ue=B'y(y+1) . . (y+r—1), (2) 
from which it follows immediately that the mean is 


= By. (3) 
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From the moment relationships the second, third, and 
fourth moments about the mean are easily found to be 


(4) 
By (5) 
3B*y(y+2). (6) 


Since the skewness statistic is Vb,=/0%, we have 
from (4) and (5) that 


(7) 


Hence, the skewness goes to zero with increasing y showing 
that the gamma distribution becomes symmetrical for 
large y; in fact, it may be shown that the distribution 
approaches normality slowly as y increases. For y>100 
it is approximately normal for climatological applications. 

The main interest in applications to climatological anal- 
ysis is not in equation (1) but in its integral which gives 
probability. This cannot be found except as an expansion 
in series or continued fractions. The integral from 0 to 
any value of the variate has been tabulated by Pearson 
[4], and this of course gives the probability that any value 
of the variate is less than the tabulated value. 

Pearson used the moments for fitting the Type III fre- 
quency curve, so the arguments of his table are w and p. 
The variate w is scaled in terms of the standard deviation 
instead of 8. Hence, u=z/(6yy) or 2/B=wy/y in our nota- 
tion. Also, his p=y—1. To use Pearson’s table we find 
y=p+l1 and multiply the w value by Vy. For certain 
purposes it may be more convenient to convert the z’s to 
u values. 


3. STATISTICAL ESTIMATORS 


The main statistical problem in the application of the 
gamma distribution to climatological data is the estima- 
tion of the parameters 8 and y from a sample record. The 
estimation problem is one of three basic problems in 
statistical analysis, and it will serve our present purpose to 
discuss the general problem briefly. 

It has long been known that there are many ways to 
estimate the parameters in a statistical equation from a 
sample of data. Two of the more common methods are 
least squares and moments. It was found by Fisher [3] 
that the various methods of estimation do not give equally 
good results in the sense that some estimates or statistics 
are more variable than others. Clearly, the best estimates 
are those which have the smallest variability. For exam- 
ple, in samples of 10 from a normal population the mean or 
expected value could be estimated by averaging the smallest 
and the largest value, or it could be estimated by averaging 
all the observations. Obviously, the latter statistic using 
all the observations should be better than that using only 
two of the observations. In fact, it has been shown that 


the variability from sample to sample for sample size 10 
as measured by the variance is twice as large when only the 
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extreme observations are averaged. The median, which 
is also an estimate of mean, has a variance about one-third 
greater than the mean for sample size 10. From this it 
may be inferred that if we use the mean range as an esti- 
mate, we in effect discard half of our data; if we use the 
median we discard one-third of it. In climatology, where 
data are scarce, the use of inefficient estimators is clearly 
to be avoided. 

Fisher [3] made a remarkable contribution to statistical 
analysis by developing a method of estimation originally 
due to Gauss which he called the method of maximum 
likelihood (M. L.). This method consists of maximizing 
what he calls the likelihood or the product of the frequency 
functions of a sample. If f(z;8,y) is any frequency fune- 
tion, the likelihood is defined to be 


M=II 


where 2, is the ith value in asample of n. To maximize this 
it is simplest to take logarithms before differentiating and 
setting to zero. This gives 


log f(21;8,7). (9) 


Differentiating partially with respect to 6 and y gives the 
M. L. differential equations 


oL_ 

10) 


Solving these gives the M. L. estimates commonly written 
as Band ¥. The M. L. estimates have certain remarkable 
advantages not always possessed by other estimates which 
will now be discussed. 

In order to assess the quality of estimators in general, 
Fisher defined three desirable properties of statistics; viz., 
consistency, efficiency, and sufficiency. These may be 
defined as follows: 

1. If an estimator or statistic is consistent it converges 
in probability to its population or parameter value. This 
may be expressed in symbols by 


P(|T,—0|< 6) >1—2; n>N. (11) 


T,, is an estimate of the parameter @ based on sample size 
n, € and are arbitrarily small quantities, and N is any 
integer. This means that 7,=0 when 7, is calculated 
from the whole population. 

2. A consistent estimate 7; is said to be more efficient 
than another consistent estimate 7, if v(7;)<v(72); 1. @ 
if the variance of 7; is less than the variance T>. An esti- 
mate is said to be efficient if it has the smallest variance 
of a class of consistent estimates. The efficiency of a0 


estimate is defined as v(7)/v(T) where 7 is M. L. estimate. 
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3. An estimate 7’ is said to be sufficient if it exhausts 
all possible information on.@ from a sample of any size. 
If T, and 7, are two different estimates of @ not func- 
tionally related, an estimate 7; of @ is sufficient if the joint 
distribution of 7, and 7, has the form 


F=A(Ti, T2|T;) (12) 


where f; is the frequency distribution of 7, and f, is the 
distribution of 7, given a sample value of 7;. Once 7; 
is known the probability of any range of values for 7; is 
the same for all 6; hence, 7; cannot give any information 
on 6 which is not already available from 7;. Sufficiency 
is the most desirable property of an estimate, and such 
estimates are said to be optimum. 

The superiority of M. L. estimates was demonstrated 
by Fisher and others when they proved that M. L. esti- 
mates are consistent and efficient and if a sufficient esti- 
mate exists, it will be given by the M. L. method. 


4. MAXIMUM LIKELIHOOD ESTIMATES 


Applying (9) to the gamma distribution equation (1) 
gives 


L=—ny log B—nlogl (y) + (y—1) log (13) 


where the summation is over the nm sample values. Dif- 
ferentiating as indicated in equations (10) we find the 
M. L. equations 


z/8—j=0 (14) 


a, 
log log r(7) log z=0. (15) 


Since x log T' (7) is the digamma function, ¥(7), we may 
write (15) in the simplified form 


log log 2=0. (16) 


Taking logarithms of (14) and substituting for log @ in 
(16) gives 


log = log log (17) 
This equation is implicit in 7 but may be solved with 
some difficulty using the Davis [5] tables of the y-func- 
tions. Masuyama and Kuroiwa [6] prepared tables of 
log ¥—y(¥) from tables of logarithms and tables of the 
digamma functions. 

We developed the application of the gamma distribu- 
tion to precipitation before Masuyama and Kuroiwa’s 
tables were available although, of course, we had also 
followed the equivalent procedure of using the Davis 
tables. To simplify the technique of fitting we developed 
al approximation to log y—y¥(7) as follows: Nérlund [7] 
shows that 
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¥(1)= log "By/(2k) + Rm (18) 


is an asymptotic expansion in which B, are the Bernoulli 
numbers, B;=1/6, B,=1/30, etc., and R,, is the remainder 
after m terms. For y>1 we may write the inequality 


Bn 1 


|Rm (19) 


For only m=1 and y=1, |R,,| <0.00833 which is less 
than 1.5 percent of the table value y (1) = —0.57722 given 
by Davis [5]. The approximation, of course, increases 
in accuracy with y. At y=2 it is within 0.1 percent of 
table value. We are not, however, interested in approxi- 
mating y but in approximating y. 

From (18) for m=1 we find 


¥(y) =log y—1/(2y) —1/(12y’). (20) 
Substituting in (17) we find 
12 (log log ) 7*—6y—1=0. 
Simplifying by letting A=log z—* = log z we have 
12 (21) 


which is a quadratic equation whose only pertinent root is 


This together with equation (14) gives the M. L. estimates 
It is only necessary to sum 
the natural logarithms of z and take the natural logarithm 
of the mean of z to provide the basic data for equations 
(14) and (21). Common logarithms may, of course, also 
be used by multiplying by the proper conversion factor. 
The error in 7 resulting from using only one term of 
equation (18) is not readily expressed in mathematical 
form; hence, we have computed the following table for 
correcting the estimate obtained from equation (22). 


(22) 


0.2 0.034 0.8 0.012 1.4 06.006 2.2 0. 003 
0.3 .029 0.9 .O11 1.5 005 2.3 002 
0.4 .025 10 .009 1.6 . 005 3.1 . 002 
0.5 021 1.1 .008 1,7 . 004 3. 2 . 001 
0.6 .017 12 .007 1.8 . 004 5.5 . 001 
0.7 .014 13 .006 1.9 . 003 5. 6 . 000 


The value of Ay is to be subtracted from the value of 7 
obtained from equation (22). 


5, SUFFICIENCY OF THE ESTIMATES 


Koopman [8] has given the necessary and sufficient 
conditions for a set of estimators to be jointly sufficient 
in the form 
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log f= (23) 


Here A, and B are functions of the parameters of the 
distribution (8 and y in this case) and X, and Y are func- 
tions of z. Taking logarithms of (1) we have 


log f= —zx/B+(y—1) log z—log T'(y)—vy log 8. — (24) 


This is of the form (23) where A;,=—1/8, A,=y—1, 
X\=2, X,=log Y=0, and B= (log! (y) +7 log B). 
This shows that # and 7 are jointly sufficient estimates 
of 6 and y. Thus, no other estimates of 6 and y can 
give more information on these parameters, and # and 7 
are optimum statistics and have, indeed, a highly desirable 
property. Fisher [3] has also shown that Z is a 100 percent 
efficient estimate of the population mean. This is im- 
portant in many applications of the gamma distribution. 


6. VARIANCE OF THE ESTIMATORS 


Fisher [3] has shown how to obtain the large sample 
variance of M. L. estimates in general and has demon- 
strated that M. L. estimates are normally distributed in 
large samples. Using matrix methods we define 


ES log (25) 


where o” is the ijth element of a matrix, the 6’s are the 
parameters, and £ is the expected value operator. If 
further we define the matrix equation 


[o. +9] [o (26) 


the variance-covariance matrix is 


(27) 


Since 6 takes on two values 6 and y, we may write the 
second partial derivatives as Z,;, L2, and Ly. where the 
subscripts refer to the parameters. Differentiating equa- 
tion (24) partially with respect to 6 and y we find 


Y 2 
(28) 
Ly= —y (Y) (29) 
Ly= In= (30) 


where y’ is the trigamma function. 
We have seen above in equation (3) that (x) 
hence 


(31) 


Since (29) and (30) do not involve the random variable z, 
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we have simply 
(7) (32) 
1 
We now may write the matrix 
a 1 
8B 
[o]= (34) 


The inverse of this after multiplying each term by 1/n as 
required by (27) is 
By’ (7) 
—B 
(Y)—1) 


(35) 


which is the variance-covariance matrix. From this it 
immediately follows that the variances and covariance are 


(6) 


(38) 


Since cov (8,7) =ryv(§)0(7), we easily find the correlation 
r between § and 7 to be 


—1 


This is negative as it is seen it must be by equation (14), 
for if the mean is fixed, 8 and 7 must vary inversely. The 
large sample variances were also given by Masuyama and 
Kuroiwa. 


7. EFFICIENCY OF THE METHOD OF MOMENTS 


Although the method of moments is widely used in 
fitting frequency distributions in climatology, it is the 
exceptional case when this method proves to be fully 
efficient in estimating climatological parameters. The 
most prominent exception is, of course, the normal dis- 
tribution where the moment estimates are jointly sufficient 
and so are identical with the M. L. estimates. As we 
have seen the M. L. method of estimation is always 
superior to moments when they give different results. 

It is of interest then to compute the efficiency of the 
moment estimates for the gamma distribution to evaluate 


fo 
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the superiority of M. L. estimates. We have seen that 
the efficiency of a moment estimate is the ratio of the 
variance of the M. L. estimate to that. of the moment 
estimate. The M. L. variances are given by equations 
(36) and (37) and the variances of the moment estimates 
may be worked out as follows: 


(40) 
Ym = 2/82. (41) 


The problem now is to find the variance of these two 
estimates. This may be done using Taylor’s formula. 
Taking the total differential of (40) with respect to Z 
and gives 


AF (42) 
z z? 
in which the coefficients of the delta increments are the 


partial derivatives and the increments of higher order are 
ignored. Squaring and taking expected values we find 


E(asy +5, E(Az E(ATAs*). (43) 


Since the expected value of the square of an error is the 
variance, and the expected value of an error product is 
the covariance, we have 


(6) +55 cov (Zs). (44) 

Now the variances ofs?, Z, and cov (7, s”) are known in large 
samples to be [9] 

v(s*) = sa") (45) 

o(Z) = (46) 

cov (Z, s?)=(n—1) (47) 


Substituting values of the moments of the gamma dis- 
tribution given by equations (3), (4), (5), and (6) for_the 
moments in (45), (46), and (47) and then for the sample 
moments, variances, and covariance in (44) we find after 
some simplification 


0(6m) (48) 
which for large n becomes 
wt 
0(Bm) (2y+3). (49) 


Performing similar operations on equation (41) we find 
for large samples 


==" (y-+1). (50) 
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If we assume 8=1, a common value for climatological 
data, we find for y=1 | 


(Bm) =5/n; 0(¥m)=4/n. (51) 
and for y=10, 
0(Bm)=2.3/n; v (¥m)=220/n. (52) 


Using the Davis tables [5] we find the M. L. variances 
from (36) and (37) for y=1 


v(§)=2.55/n; v (4)=1.55/n (53) 
and for y=10, 
v(8) =2.04/n; v (7) =193.6/n. (54) 


According to our definition of efficiency we must take the 
ratio of the variance of the M. L. estimate to the moment 
estimate. This gives from equation (51) and (53) for 
y=1, 


Eff (8,,)=51 percent 
Eff (ym)=39 percent 


Eff (8,.)=89 percent 
Eff (7¥m)=88 percent 


and for y=10 


We see that for y < 10 the method of moments produces 
unacceptable estimates for both 6 and y. For y near 1 
the moment estimates use only 50 percent of the informa- 
tion in the sample for estimating 8 and only 40 percent 
for y. Thus, for 6 the M. L. estimator would do as well 
with half the length of record as the moment estimate 
and for y with two-fifths the record length. Hence, use 
of the moment estimates in effect results in discarding 
half the record in estimating 8 and three-fifths of the record 
in estimating y. For y=10 the efficiencies of 8, and ym 
both approach satisfactory values. In view of the diffi- 
culty in obtaining homogeneous climatological data series, 
it seems at least a questionable procedure to employ in- 
efficient estimators which do not make the best use of the 
available climatological samples. 
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A GRAPHICAL TECHNIQUE FOR DETERMINING EVAPOTRANSPIRATION 
BY THE THORNTHWAITE METHOD 
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U. S. Weather Bureau, Washington, D. C. 


and 


A. VAUGHN HAVENS 
Department of Meteorology, Rutgers University, New Brunswick, N. J. 
[Manuscript received December 4, 1957; revised March 6, 1958] 


ABSTRACT 


Generalized diagrams for the determination of evapotranspiration and potential evapotranspiration in accordance 
with the empirical equations of Thornthwaite are presented, 


1. INTRODUCTION 


The loss of water from the earth to the atmosphere by 
transpiration from vegetation and by direct evaporation 
constitutes an important part of the water balance 
problem. However, direct measurement of these factors 
has proved to be extremely difficult, and this inherent 
difficulty has led to the development of a number of 
formulas designed to estimate water loss directly from 
meteorological data. 

It is not the purpose of this article to review or evaluate 
the various methods which have been developed, but it 
would appear that some mention of the various approaches 
to the problem is in order. A detailed and excellent 
review has been published by van Wijk and de Vries [1]; 
and an enlightening analysis of the problems involved 
was presented by Penman [2]. - 

The list of references following this article includes some 
of the principal publications devoted to the problem of 
estimating soil moisture loss directly from meteorological 
data and will suffice to indicate the nature of the various 
approaches used. They fall into three groups: (1) Those 
involving the flux of water vapor; (2) those utilizing the 
heat balance of the evaporating or transpiring surface; 
and (3) those which use an empirically determined rela- 
tionship between evapotranspiration and one or more of 
the meteorological factors involved. 

None of these methods provides a completely adequate 
solution to the problems of evaporimetry because none is 
free from assumptions, arbitrary constants, or technical 
difficulties of observation and measurement. In spite of 
the shortcomings, a number of workers have surmised that 
these methods enable the climatologist to estimate total 
évapotranspiration from a sizable field more accurately 
than a soil scientist can measure it. 


The Thornthwaite method [3] was developed from rain- 
fall and runoff data for several drainage basins. The 
result is basically an empirical relationship between 
potential evapotranspiration and mean air temperature. 
In spite of the inherent simplicity and obvious limitations 
of the method, it does surprisingly well. It is not neces- 
sarily the most accurate method, nor does it have the 
soundest theoretical basis. On the contrary, these 
distinctions probably belong to one of the vapor flux or 
heat balance methods. Among the more obvious short- 
comings of Thornthwaite’s empirical relationship is the 
inherent assumption that a high correlation exists between 
mean temperature and some of the other pertinent param- 
eters such as radiation, atmospheric moisture, and wind. 
While such limitations may be relatively unimportant 
under certain circumstances, they are at times of the 
utmost importance; and Thornthwaite and Mather [4], 
recognizing that solar radiation and atmospheric turbu- 
lence are the important factors in natural evaporation, 
have stated that the problem of developing a formula for 
potential evapotranspiration remains unsolved. 

Although ease of application is not a suitable criterion 
of adequacy, it is often a primary consideration for use. 
The vapor flux and heat balance methods require meteoro- 
logical data which are either not observed or are observed 
only at a few widely separated points. On the other 
hand, Thornthwaite’s empirical formula can be used for 
any location at which daily maximum and minimum tem- 
peratures are recorded. It is this simple universal appli- 
cability rather than any claim to outstanding accuracy 
which has led to the widespread use of this method. 

It is not the purpose of this paper to make a case either 
for or against this or any other method. Rather, recog- 
nizing the widespread interest in computing and testing 


36, 
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results of Thornthwaite’s method, it is the purpose to 
present a graphical solution which will make such com- 
putations easier. 


2. GRAPHICAL SOLUTION OF THORNTHWAITE 
METHOD 


Thornthwaite’s basic formula in the form for computing 
monthly potential evapotranspiration is 


e=1.6 (107/7)* 


where e=monthly potential evapotranspiration (cm.) 

T=monthly mean temperature (° C.) 

I=a heat index which is a constant for a given 
location and is the sum of 12 monthly index 
values 7, where 7 is a function of the monthly 
normal temperatures. (See table 1.) 

a=an empirically determined exponent which is a 
function of J, '—7.71 
+1.79X10-? [+-0.49. 


Arithmetic solution of the equation becomes an extremely 
laborious procedure primarily because of the complexity 
of the exponenta. Therefore, extensive use of the method 
for many locations over a long period of time requires the 
use of a tabular or graphical solution to the equation. 
Thornthwaite and Mather [4] have published such tables 
for only one location, Seabrook, N. J. Their publication 
includes a diagram for the graphical determination of un- 
adjusted potential evapotranspiration in metric units. 
However, the day-length correction, depending on date 
and latitude, appears in a separate table. It was thought 
that a generalized and complete graphical solution in 
English units would be useful. Such asolution is presented 
here. It includes the adjustment for day length and pro- 
vides a means whereby tables or graphs for use at any 
particular place can be readily prepared. In addition 
diagrams are provided for the conversion of weekly rates 

to daily rates of potential evapotranspiration and for a 
correction due to soil dryness. 

Weekly rates of unadjusted potential evapotranspiration 
may be obtained from figure 1, using only the mean tem- 
perature for the period and the heat index (J) for the 
location concerned. The graph also provides the adjust- 
ment for day-length required to convert unadjusted poten- 
tial evapotranspiration to adjusted potential evapotrans- 
piration. 

Figure 2 is a simple graph for converting weekly rates of 
adjusted potential evapotranspiration (PE) to daily rates. 
For certain purposes, daily rates may be desirable, 
although some evidence has been obtained which indicates 
that weekly rates of PE correlate better with measured 
evapotranspiration than do daily rates. A monthly rate 
can of course be obtained by simply multiplying the 
daily rate by the number of days in the month. 

There is considerable evidence that during long dry 
periods the rate of evapotranspiration decreases as the 
soil dries. Thornthwaite and Mather [4] suggested one 


method of correcting for this effect. It is shown graphi- 
cally in figure 3. 

Listed below is the step-by-step procedure to be fol- 
lowed in the use of these graphs. 


A. To obtain potential evapotranspiration, PE: 
1. Using monthly normal or long-term mean temper- 
atures (F.°) for the station or area concerned, 
obtain from table 1 the monthly heat index, i, 
corresponding to the normal temperature for 
each of the 12 months. 
2. Add these 12 7-values to obtain I, the heat index. 


3. Enter figure 1 at the bottom with the mean temper- 
ature (F.°) for the period concerned. (The 
length of the period is immaterial at this point; 
it may be of any length from 1 day to 1 month.) 

. Follow vertically up the appropriate mean temper- 
ature line until it intersects the horizontal heat 
index line equal to the heat index determined in 
step A-2 above. At this intersection determine 
unadjusted PE (inches) from the slanting lines. 
This is the potential moisture loss for a 7-day 
period, each day having a length of 12 hours. 

5. Enter the bottom (or top) of the upper-right por- 

tion of figure 1 with the unadjusted PE value 
found in step A-4. This determines the appro- 
priate vertical line to follow in this portion of the 


figure. 

6. Enter the upper-left portion of figure 1 according 
to the north latitude of the station or area under 
consideration and proceed vertically downward 
to the intersection with the “month” line corre- 
sponding to the month under consideration. 
This intersection determines the appropriate 
horizontal line to follow into the upper-right 
portion of figure 1. (The climatological week 
numbers and the day-length correction are dis- 
cussed in Appendix 1.) 

7. At the intersection of the vertical line determined 
in step A-5 and the horizontal line determined in 
step A-6 read the final PE in inches per week 
from the curved lines. 

8. If daily values are required, the daily mean temper- 
ature would be used in step A-3. The adjusted 
weekly PE rate obtained from step A-7 would 
then be converted to daily PE by dividing by 7 
or by using figure 2. 

9. If monthly PE is desired, monthly mean temper- 
ature would be used in step A-3 proceeding 
through step A-8 to obtain the mean daily rate 
for the month. This may be multiplied by the 
number of days in the month to obtain total PE 
for the month. 

B. To correct PE to account for soil moisture depletion: 

1. From actual measurement or from water balance 
tabulations, obtain the soil moisture content. 
This should be expressed as a percentage of the 
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Figure 1.—Nomogram for computation of potential evapotranspiration (PE). 
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Tasie 1.—Monthly values of i, according to monthly normal or long-term mean temperature (to tenths). To obtain I for use in figure 1, add 
the i-values obtained here for each of the 12 months. 

Monthly Monthly 
normal 0 1 2 3 4 5 6 x 8 9 normal 0 Jl 2 3 4 5 6 a 38 9 
temperature perat 
F°) (F°) 
0 0 0 0 0.01 | 0.01} 0.01] 0.02] 0.02] 0.03 7.49) 7.52| 7.55] 7.58| 7.61| 7.64] 7.68] 7.72] 7.75] 727 
0.08 | 0.04] 0.04| 0.05] .06| .06| .07] .08| 09 7.82) 7.85| 7.88| 7.91| 7.94] 7.98] 802| 810) gy 
.31] .32] .34] .35] 1.36] 1.37] .38| .39 
42} 143] 144] 145] 146] 147] 148] 249] |] 885 | 888| 891| 896/ 9.00) 9.04] on 
.66| .68 9.20] 9.23| 9.28] 932] 9.36] 939] 942] 945] 949) om 
9.93 | 9.97 | 10.01 | 10.04 | 10.08 | 10.11 | 10.15 | 10.18 | 10.22 | 10,% 
.85| .86| .91] .94| .95] .97| .98 30 | 10.35 | 10.38 | 10.41 | 10.44 | 10.48 | 10.52 | 10.56 | 10.60 | 10,64 
1.00] 1.01] 1.03] 1.04] 1.06] 109] 1.10] 1.12] 114] 1.16 68 | 10.71 | 10.75 | 10.78 | 10.82 | 10.85 | 10.89 | 10.94 | 10.98] 1m 
117| 1.19] 120] 1.22] 1.25] 127] 1.29] 1.30] 1.33 | 11.09 | 11.12 | 11.16 | 11.19 | 11.23 | 11.28 | 11.32 | 11.37| 
1.35 | 1.37] 1.39| 1.42] 143] 145] 147] 149] 150] 1.52 | 11.47 | 11.50 | 11.54 | 11.58 | 11.62 | 11.66 | 11.70 | 11.74] 
1.54| 1.56| 1.59] 1.61] 1.63| 1.64] 1.68| 1.70] 1.72 | 11.86 | 11.90 | 11.94 | 11.98 | 12.02 | 12.06 | 12.10 | 12.14] 1218 
1.74] 1.77] 1.79| 181] 1.83] 1.85] 1.87] 189] 191] 193 22 | 12.26 | 12.30 | 12.34 | 12.37 | 12.41 | 12.45 | 12.49 | 12.53] 12.57 
196 | 1.98} 200] 202] 204] 210] 212] 215 
217| 2.19| 2.21| 223] 225] 227/ 229) 235| 237 12.65 | 12.69 | 12.73 | 12.77 | 12.81 | 12.85 | 14.89 | 12.93] 129 
2.30) 244| 246] 248] 250| 252] 256| 258] 2.60 13.05 | 13.09 | 13.13 | 13.17 | 13.21 | 13.25 | 13.29 | 13.33] 1397 
2.62| 264] 2.66) 269] 271) 273| 279| 281] 2.83 13.45 | 13.49 | 13.53 | 13.57 | 13.61 | 13.65 | 13.69 | 13.73 | 13.7 
13.85 | 13.89 | 13.93 | 13.97 | 14.01 | 14.05 | 14.09 | 14.14] i419 
2.86} 288| 292] 294] 299] 3.03/ 305] 3.08 14.28 | 14.32 | 14.36 | 14.40 | 14 44 | 14.48 | 14.52 | 14.57] 1461 
3.10} 3.12] 314| 3.17) 320] 323) 325| 327] 330) 3.32 14.69 | 14.73 | 14.77 | 14.81 | 14.85 | 14.90 | 14.95 | 14.99] 15.6 
3.34] 3.36| 3.30) 342] 3.48] 353| 355] 3.58 45.11 | 15.15 | 15.19 | 15.23 | 15.28 | 15.33 | 15.38 | 15.42] 15.4 
3.60} 3.62] 3.67| 3.60] 3.72) 3.74| 3.76| 3.78] 3.83 15. 54 | 15.58 | 15.62 | 15.66 | 15.70 | 15.75 | 15.79 | 15.84 | 15.8 
3.86 | 3.90| 3.93) 3.96] 3.98] 400] 402] 405| 407] 410 15.96 | 16.01 | 16.06 | 16.10 | 16.15 | 16.19 | 16.23 | 16.27] 16.3 
414] 417| 422| 425] 427| 430] 432] 439 16.39 | 16.44 | 16.49 | 16.54 | 16.58 | 16.62 | 16.66 | 16.70 | 16.7% 
4.42| 445| 447] 4.50] 452] 455| 457| 460| 464] 467 
Tisedekenth 4.70} 4.72| 475| 4.79] 4.81] 483| 486] 490] 493| 496 16.83 | 16.88 | 16.93 | 16.97 | 17.01 | 17.05 | 17.09 | 17.13 | 17.18 
4.98} 5.01| 5.04] 5.07] 5.09] 512] 5.16] 5.19| 5.22] 6.25 17. 28 | 17.33 | 17.37 | 17.41 | 17.45 | 17.49 | 17.53 | 17.58 | 17.8 
5.28 | 6.30] 5.33| 5.35] 5.39] 5.49] 5.52] 5.55 17.73 | 17.77 | 17.81 | 17.85 | 17.89 | 17.93 | 17.98 | 18.03 | 1808 
18.17 | 18.22 | 18 26 | 18.30 | 18.34 | 18.39 | 18.44 | 18.49] 18 
5.57 | 6.60] 5.63| 5.66] 5.70| 5.73| 5.76| 5.79] 5.82| 5.84 18. 62 | 18.66 | 18.70 | 18.75 | 18.80 | 18.85 | 18.90 | 18.95 | 18% 
5.87 | 5.90| 5.93| 5.98] 6.01| 6.04) 6.07] 610) 612] 6.15 19.07 | 19.12 | 19.16 | 19.21 | 19.26 | 19.31 | 19.36 | 19.41 | 19.45 
6.17| 621| 626| 629] 635 | 638| 641| 644] 6.47 19.53 | 19.58 | 19.63 | 19.68 | 19.73 | 19.78 | 19.82 | 19.86 | 
6.50 | 6.54] 6.57| 6.60] 6.63 | 6.66| 6.69| 6.72] 6.75| 6.78 20.00 | 20.05 | 20.10 | 20.15 | 20.20 | 20.24 | 20.28 | 20.32) 2.97 
6.84| 6.81/ 6.90} 692] 695 | 698| 7.01] 7.04] 7.08] 7.12 20.47 | 20.52 | 20.57 | 20.62 | 20.66 | 20.70 | 20.75 | 20.80 | 20.8% 
7.15 | 7.19| 7.22] 7.26] 7.28] 7.31] 7.34] 7.38] 7.42] 7.45 20.95 | 21.00 | 21.04 | 21.09 | 21.13 | 21.17 | 21.21 | 21.25 | 
total water-holding capacity in the root zone of of field capacity. (Since the heat index (J) is given, steps 


the crop when soil moisture is at field capacity. 

2. Enter figure 3 from the bottom with the daily PE 
value obtained from step A-8. Proceed upward 
along this line to the appropriate diagonal line 
representing the percentage value obtained in 
step B-1.. Read the corrected evapotranspiration 
(E) along the diagonal scale at the left. 


3. EXAMPLE 


Find the evapotranspiration (E) for July 15 (climato- 
logical week No. 20) at a station located at 41° N. and 
having a heat index (J) of 52. Assume a daily mean tem- 
perature of 70° F. and a soil moisture content of 70 percent 


A-1 and A-2 are not necessary.) 

1. Proceed through step A-4, using the given mean 
temperature of 70° F. and J value of 52. This 
gives a weekly unadjusted PE rate of 0.92 inch. 

2. Using the latitude of 41° N. and the date of July 15 
proceed through steps A-5 to A-7 and determine 
the adjusted weekly PE rate of 1.13 inches. 

3. To convert this weekly rate to a daily value pro- 
ceed through step A-8 using figure 2. The final 
value is about 0.16 inch per day. 

4. Finally, to correct this PE rate of 0.16 inch per day 
for the effect of soil moisture depletion, using the 
70 percent of capacity as given, proceed through 
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Ficure 2.—Nomogram for conversion of weekly PE to daily PE. 
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Figure 3.—Nomogram for adjustment of PE for scil moisture 
depletion. 


step B-2 using figure 3. A daily rate of about 
0.11 inch is determined, which is now an esti- 
mate of actual evapotranspiration rather than 
potential. 


4. CONCLUSIONS 


This graphical method is most useful for work involving 
a wide range of climatic conditions (i. e., a wide range of J 
values.) For computations covering a long period of time 
at a single station or a group of stations with nearly equal 
I and latitude values, it will usually save time to use 
tables or graphs which apply only to those particular 
stations. Such tables or graphs may easily be prepared 
using the diagrams presented here over a suitable range of 
temperature and time. For example, at a given place, 
latitude and heat index are constants and PE depends 
only on temperature and time of year; so a table of PE 
as a function of these two variables could be prepared from 
figure 1 without difficulty. (See Appendix 1.) 


APPENDIX 1 


The upper-left portion of figure 1 concerns the correc- 
tion for latitude and time of year. The unadjusted PE 
obtained in the lower portion of figure 1 must be multi- 
plied by this day-length correction in order to obtain 
adjusted rate of weekly PE. The curved lines labeled 
by months represent mid-month, and for most purposes 
interpolation between these lines according to date will 
suffice. However, for a more accurate determination 
one may use climatological week number (March 1-7== 
week No. 1, 8-14—No. 2, etc.) rather than month. These 
lines are not reproduced, but the week numbers are shown 
for the middle day of the week at 30° and 50° N. latitude. 
With an appropriate drafting curve one can construct the 
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Ficure 4.—Relation of heat index, I, to rormal annual temperature. 


lines for week number. For example, the curve for week 
No. 10 (May 3-9) would have approximately the same 
shape as the May line, whereas week No. 24 (Aug. 16-22) 
would approximate the August curvature. 

It is not necessary that the numerical value of the day- 
length correction be determined; the adjustment is ac- 
complished graphically when the appropriate horizontal 
line of day-length correction is followed to the right until 
it intersects the appropriate vertical line of unadjusted 
PE and the final PE is read from the curved lines. How- 
ever, in order to facilitate the work of those who may wish 
to prepare tables for a particular station, the scale for day- 
length correction has been included. 

The preparation of tables for any given location is a 
simple matter. A table of unadjusted PE as a function 
of temperature can be prepared from the lower portion 
of figure 1 along the appropriate heat index line. The 
day-length correction table can be prepared from the 
upper-left portion of figure 1 by tabulating the value of 
the correction for each month (or week) along the appro- 
priate latitude line. The final PE for any particular week 
isthen the product of the appropriate two values from the 
tables. Another convenient device for use at a particular 
station is a set of tables or graphs of final PE as a function 
of temperature for each of the 12 months. If these tables 
are made for monthly PE, an adjustment for month length 
must be made as indicated in step A-9. 

Figure 1 as drawn applies to the Northern Hemisphere. 
It could be adapted for use in the Southern Hemisphere 
by the appropriate relabeling of the month lines in the 
upper-left portion. 


APPENDIX 2 


Inasmuch as monthly heat index, i, is a function of 
temperature only [3] 


i=(T/5)=™, 


it seemed reasonable to investigate the relationship be- 
tween heat index, J, as determined from table 1, and 
normal annual temperature. This was done for the 43 
stations shown in the accompanying table 2. The only 
criterion for the selection of these stations was that a 
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TABLE 2.—Normal annual temperature and heat index value, I. (Arranged in order of increasing temperature.) 
ann 4 ann ex I, 
Station tempera- from table 1 Station pera- from table} 
ture (°F.) ture (°F.) 

53.5 55 || Ponce Santa Isobel, West 76.7 136 
54.7 58 || Swan Island, West Indies. 80.7 1% 
wide and uninterrupted range of normal annual tempera- 2. H. L. Penman, “Evaporation: An Introductory Survey,” 


tures should be represented. As a consequence a number 
of climatic types are represented but not equally so. 
Figure 4 shows heat index plotted as a function of nor- 
mal annual temperature for these 43 points. The curve 
was drawn by eye. The relationship is surprisingly good 
and suggests that heat index, J, can be quickly estimated 
from normal annual temperature with a maximum error 
of + 4 or 5. This amount of error in the heat index 
would introduce an error in PE of no more than .03 or 
.04 inch per week, which is very likely well within the 
limits of accuracy of the entire empirical procedure. 
Of course, this result suggests that the heat index scale 
in figure 1 could be replaced by a scale of normal annual 
temperature without significant loss of accuracy. 
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RADIATIONAL INVERSIONS AND SURFACE TEMPERATURE CHANGES 


SIGMUND FRITZ 
U. S. Weather Bureau, Washington, D. C. 
[Manuscript received March 12, 1958; revised April 24, 1958] 


ABSTRACT 


When air stagnates in dark polar regions, an inversion forms. The change of the inversion magnitude is studied 
under the assumption that the snow surface radiates about as much energy as it receives from the atmosphere. It 
turns out that the inversion magnitude may either decrease or increase as the surface temperature falls, depending 
on the rate of change of atmospheric “emissivity” with air temperature, 


1, INTRODUCTION 


Wexler [1, 2] has studied the formation of inversions 
during the long polar night. An examination of his 
papers shows that under some conditions the inversion 
magnitude, AT, decreases, and for other conditions AT’ 
increases as the surface temperature, J, decreases. The 
inversion magnitude is given by 


AT=T,—T, (1) 


where 7, is the temperature of an isothermal layer of air 
which lies above a colder surface. 

Actually, Wexler studied the following phenomenon. 
He assumed that during the long polar night a relatively 
unstable maritime airmass stagnated over a snow surface 
with initial surface temperature 273° K. He then tried 
to compute the future course of the air temperature, 7, 
and surface temperature, 7,, assuming that only radiation 
was acting. The model he adopted stipulated that 7, 
falls very quickly to a quasi-equilibrium value determired 
by the downcoming radiation from the air. Thereafter, 
T, decreases slowly. The air temperature then also falls 
slowly, forming an isothermal layer with temperature, 
T., above the ground; above the isothermal layer, the 
original maritime air temperature prevails. In his 1936 
paper, with the aid of water vapor absorption coefficients 
measured by Weber and Randall [3], Wexler [1] found 
that AT decreased as T, decreased, after 7, had reached 
the quasi-equilibrium value. Later, in 1941, Wexler [2] 
found that if he used the newly developed Elsasser diagram 
for saturated air AT increased as 7, decreased. Whether 
or not AT increases or decreases depends critically on the 
rate of change of the atmospheric radiation with air 
temperature. 

It is the purpose of this paper to inquire under what 
conditions AJ’ would increase with decreasing 7’, and 
under what conditions AJ’ would decrease with decreasing 
T,, in an individual mass of air cooling by radiation only. 


2. EMISSIVITY AND RADIATION 


From the assumed state of quasi-radiative equilibrium 
between the snow surface and the overlying air, we get 


R=cT,'=koT,' (2) 
T,.='“T, (3) 


where & is an overall “emissivity” of the atmosphere, and 
we assume that the main part of the downward atmos- 
pheric radiation comes form the isothermal layer; o is 
the Stefan-Boltzmann constant. 

Since AT=T,—T;,, by substituting for 7, from equa- 
tion (3), we get 


AT=T,(1—k') (4) 


The problem now is to find under what conditions of radia- 
tion, R, and therefore of atmospheric emissivity, k, AT’ 
increases with decreasing 7, and when AT’ decreases with 
decreasing 7T,. We can get the limiting condition by 
making AT’ constant. 

Thus equation (4) becomes 


4.T 
qe (5) 


where A, denotes that AT is constant with variations in 
T, or T,.* Equation (5) shows the relation between the 
““emissivity’’, k, and isothermal air temperature, 7,, which 
is required in order for the inversion magnitude, A7’, to 
remain constant as the surface temperature, 7’, drops. 
This relation of & vs. T, is shown in figure 1 for a range of 
k and T, which is found in Wexler’s studies. Figure 1 
shows that the emissivity of the isothermal air may de- 
crease with decreasing air temperature even though the 


*Under the conditions that 47'=constant with 7’o, 


dT,” oT. 
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Ficure 1.—Variations of atmospheric “emissivity,” k, with inversion 
temperature, 7, when the inversion magnitude, A.7, does not 
change with surface temperature, T7,. 


inversion magnitude is not changing as surface temperature 
falls. However, if the slope of the k vs. T, lines is steeper 
than those in figure 1 (that is, if k decreases faster with 
decreasing 7), AJ’ will increase as 7’, decreases; if the 
slope is less steep, AT’ will decrease. 

The conditions on the atmospheric radiation R, for 
AT=constant, can be readily found from equation (2). 
As background, however, we should look at 2 vs. T curves 
from Wexler’s 1941 paper [2]. In figure 2, curves a, c, 
and d are taken from that paper. Curve a represents the 
black body radiation of a snow surface, plotted against 
T,; curves c and d represent the downward atmospheric 
radiation computed in various ways and plotted against 
T,. In curve ec, 7, represents “the highest mean tem- 
perature of any layer of air containing 1 mm. of precipi- 
table water and normal CO, content.” In the computa- 
tion of R for curve c the absorption coefficients of Weber 
and Randall [3] were used. Curve d was “. . . computed 
by use of the Elsasser diagram in the following way: the 
radiation coming from a saturated atmosphere in convec- 
tive equilibrium with an ocean surface of 0° C. is desig- 


refer to the radiation coming from atmospheres whose 
lower portions are isothermal at temperatures correspond- 
ing to the abscissae of the curve and whose upper portions 
follow the original convective equilibrium curve.” The 
values of 7, and 7’, tabulated by Wexler [2] show that 


> 0 for curve c and less than zero for curve d.** 


This result is to be expected. From the definition of 
A.T, by substituting for 7, in equation (2), we get 


R=o(T,—A.T)* 


The values of RP vs. 7, for several values of A,7 are plotted 
in figure 2. These curves, of course, apply only when AT 
is constant with changes in 7, or 7,. If any lines in 
figure (2) converge toward curve a, as TJ, decreases, more 


rapidly than the A,7' curves, then a >0; if such lines 


diverge from curve a more rapidly than the A,7’ curves 


then ue <0. Thus, curves c and d represent the type 


of variation of oar to be expected; that is, for curve ¢ 


the inversion decreases with decreasing surface temper- 
ature; the opposite is true for curve d. It is easy to 
verify from figure 2, that in curve c, k is very nearly 
constant. But for curve d, k decreases more rapidly with 
decreasing 7, than indicated by the curves in figure 1. 


**Wexler’s table shows this for T.<—40° C. However his curve d shows it for al] 
values of Ts. 
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Thus, we might expect in our simplified model, the 

condition for 227 <0 is that the “emissivity” of the 


atmosphere should decrease more rapidily with decreasing 
air temperature than the curves in figure 1. 


3. DISCUSSION 


Although curves c and d represent only computation 
with special assumptions and conditions, it is not hard to 
visualize that in actual atmospheres the emissivity may 
vary with J, in the same way as implied in curve c or in 
curve d. For example, evaporation and condensation 
may do this. Therefore, it is not valid to generalize that 
the inversion magnitude either increases or decreases with 
decreasing surface temperature after quasi-equilibrium has 
been reached. What happens in any particular case 
depends on the radiation properties of the atmosphere 
and how they change with temperature. 

Moreover comparison of 7, with AT with the aid of 
many different radio-soundings for many different meteor- 
ological conditions and airmasses, does not shed any light 
on the direction of change of AT with 7, in a single airmass. 
Nevertheless Belmont [4] has concluded from his own 
findings that his results conflict with Wexler’s theoretical 
results. Belmont compiled measurements of inversions 
and found that, on the average, the magnitude of “radia- 
tion” inversions increases as T, decreases, especially for 
low values of 7. Since radiosonde data were used, the 
results do not refer to the history of a single mass of air 
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during the time when the air and surface cool. The 
results compare, instead, the temperatures of many 
different masses of air with the simultaneous surface 
temperatures. Therefore, Belmont’s results are, in any 
case, not comparable with Wexler’s results since each was 
dealing with a different problem. But it should never- 
theless be pointed out that Wexler’s results include both 
the increase and decrease of AZ with decreasing 7,, 
depending on the rate of change of atmospheric emissivity 
with air temperature. 
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1. INTRODUCTORY REMARKS 


Not since November 1957 has there been a monthly 
mean trough in the United States west of the Mississippi 
River, although in March 1958 a trough of moderate 
intensity was located along the west coast (fig. 4 of [1]). 
During April (fig. 1) that trough moved inland, perhaps 
as a result of the downstream dispersion of vorticity 
following the generation of increased westerlies and a new 
trough in the Pacific at middle latitudes. The increase 
in zonal winds in the Pacific followed the development of 
an Asiatic coastal trough in March, a feature conspicu- 
ously absent during the months of January and February 
1958 [2, 3]. 

In April, on the average, there were remnants of 
blocking over northeastern Canada, but the complete 
collapse of that blocking near the end of the month and 
its replacement by a mean trough permitted a northerly 
displacement of mid-latitude westerlies. Retreat of 
Canadian blocking toward the west and entry of the west 
coast trough into the United States were responsible for 
considerable variability in weather and circulation pat- 
terns, not only from week to week but also on a bi-weekly 
basis. Many daily, monthly, and seasonal records of 
temperature and precipitation were broken, some of 
which are listed in table 1. 


2. MONTHLY MEAN CIRCULATION AND WEATHER 
INDEX 


The monthly mean circulation for April, shown in 
figure 1, represents a transition stage over North America 
from a low index state to one of a higher index. Of 
course, that is not apparent solely from the 700-mb. 
contours. One can see that there was a general zonal 
flow in the temperate westerlies of the Western Hemis- 
phere, with the exception of the interruptions in the Gulf 
of Alaska and northwestern Canada. However, the 
departure from normal isopleths (hereafter called DN), 
shown as dotted lines in figure 1, give a more graphic 
representation of the state of the zonal westerlies. In 
this instance there was actually an easterly flow relative 
to normal north of 45° N. over most of North America. 
Thus, in spite of the apparently well-formed westerlies, 
there were some low index features, especially blocking, 
which persisted long enough to retain their identities and 
affect the month-long time-averaged mean. 


1 See Charts I-XVIT following p. 148 for analyzed climatological data for the month, 


TaBLE 1.—April monthly weather records broken in 1958 


Station Record 

Precipitation 
aa rane Greatest snowfall in 80 years, 
Wettest of record. 

Temperature 


The computed zonal index for temperate westerlies 
(35° N.—55° N. in the Western Hemisphere) for April 
was 7.9 m. p. s., only 0.4 m. p. s. below normal. But 
that represents a considerable recovery over the previous 
month whose index was 5.3 m. p. s., some 3.8 m. p. 8. 
below normal. The rise of index in April is even more 
significant if one considers that the index normally de- 
creases by 0.8 m. p. s. from March to April. 


NORTHERN HEMISPHERE CIRCULATION 


While our interest is centered principally on North 
America, it is worthwhile to examine the salient features 
of the general circulation of the whole Northern Hemis- 
phere. In figure 1 there are three major troughs at high 
latitudes and_ six significant middle-to-low-latitude 
troughs. That wave number was supported to some 
extent by the existence of high-latitude blocking, which 
can be seen as positive DN centers over Alaska, northern 
Quebec, and the Norwegian Sea. 

The outstanding feature of the DN pattern was the 
overall weakness of the gradient, except for the area from 
Scandinavia to northwestern Siberia, and another just 
east of Japan. In the latter region the abnormal strength 
of the subtropical ridge, which was longitudinally in phase 
with a deep polar trough, accounted for the fairly strong 
westerlies sweeping across the north-central Pacific. The 
DN gradient over the United States was especially weak 
and is not indicative of the variability within the month. 

The subtropical ridging in the eastern Pacific was 
accompanied by a difluent area along the United States 
west coast, to the rear of the progressive trough in the 
central United States. The difluence was partially re 
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Ficure 1.—Mean 700-mb. contours (solid) and height departures from normal (dotted) (both in tens of feet) for April 1958. High-latitude 
blocking over North America is defined by the positive height departures in the north and negative height departures to the south. 


sponsible for lower than normal 700-mb. westerly wind 
speeds across the northern United States. The slight 
tidge in the eastern United States and the residual block 
in Quebec also contributed to the weaker than normal 
westerlies. 

Most of North America was under the influence of 
tyclonic flow at the 700-mb. level as can be seen in figure 1, 
with Canada having a flow much more cyclonic than that 
in the United States. The short wave-spacing between 
the two United States troughs along 45° N. might have 
been incompatible had it not been supported by the block 
in Quebec. 


468423—58——-2 


The Atlantic was characterized by the normally sharp 
slope (vertically) of the Icelandic Low [4] and the presence 
of a strong ridge in the eastern portion. The DN field 
shows that there must have been some tendency toward 
bridging between the Canadian block and the intensi- 
fying Atlantic ridge. The European-African trough is 
normally sheared [4] with the middle-to-low-latitude 
segment west of Spain and Africa. Its position far in- 
land this month, however, was favored by the building of 
the Atlantic ridge and the Scandinavian block. 

To a great extent Asia was dominated by the deepening 
Low in north-central Siberia and its rather extensive DN 
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Ficure 2.—Mean 700-mb. isotachs for April 1958 in meters per 
second. Solid arrows represent the mean 700-mb. jet stream 
which was located far to the south. Centers of fast and slow 
wind speed are located by F and §. 


field. The motion of that Low, some 15° of longitude 
eastward from the previous month, may have repelled 
the next Low and trough downstream, which moved a 
like distance from March to April. 

Most of the features discussed above are further illus- 
trated in figure 2, the monthly mean 700-mb. isotachs. 
Note especially the split of the jet stream into high- and 
low-latitude branches in those areas where blocking 
predominated. 


UNITED STATES WEATHER 


The departure from normal of average temperature 
in the United States for April, shown in Chart I-B, 
reveals a broad band of below normal temperatures from 
the Pacific Northwest to the Rio Grande and Tennessee 
Valleys. Except for Nevada and the Texas Panhandle, 
the magnitude of the departures was quite small, which is 
consistent with the rather flat negative DN field shown in 
figure 1. The normal to above normal temperatures 
from eastern Montana to the Middle Atlantic States 
resulted from the persistent Canadian blocking (and 
resultant easterly DN flow) and above normal 700-mb. 
heights which extended from the Northern Plains States 
eastward. A further factor to consider is that the tracks 
of polar anticyclones (Chart IX) were principally north 
of the Canadian border. 

The above normal temperatures observed in the Cali- 
fornia coastal belt cannot be explained in terms of DN 
flow. In fact, an onshore DN and contour flow would 
usually produce cooler than normal temperatures. How- 
ever, cold air inland and its suppressing effect on the sea 
breeze usually does produce some warming along the coast. 
But since warmer than normal temperatures have been 
extremely persistent in this area, there may be other 
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Ficure 3.—Change of 700-mb. mean height departures from nor. 
mal from March to April 1958. Labels of centers are in tens of 
feet. Note increased westerly flow in most of the Westem 
Hemisphere at middle latitudes. 


factors contributing to this complex subject, such as the 
warmer than normal sea surface temperatures reported 
off California for the past several months. 

April was a rather wet month in many areas of the 
United States. (See Charts II and III-B.) At least 
half of the country received greater than normal amounts 
of precipitation, with many areas exceeding more than one 
and one-half times normal. Several stations reported 
record April rains (table 1), and heavy snow fell in parts 
of the Sierra Nevada mountains, Utah, New Mexico, and 
Montana. Only one station reported a record dry month. 

In terms of the circulation, the heavy precipitation 
west of the Continental Divide fell with the advance inland 
of the west coast trough. The heavy precipitation from 
Texas eastward and northeastward closely paralleled 


‘the axis of low pressure seen on the mean sea level chart 


for April (Chart XI). Also associated with that band of 
heavy rainfall were the displacement of the 700-mb. 
jet stream far to the south (fig. 2) and the location of a 
secondary maximum of cyclone frequencies (shown only 
as tracks in Chart X) in the Gulf of Mexico, where cyclone 
frequency is normally low during April [5]. 

The area of below normal precipitation from souther 
Utah to Missouri is difficult to explain. Most qualitative 
considerations indicate the probability of heavy precipita- 
tion in that area, but a survey of daily sea level charts 
shows that there was considerable northerly flow which 
would have been discouraging to heavy precipitation. 
The mean southerly flow at sea level on Chart XI appeals 
to have been too weak to permit the establishment of 8 
pronounced moist tongue of Gulf air that far to the north. 
This and other features mentioned above may have mort 
meaning after a discussion of the changes in circulation 
which took place from March to April. 
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3. MARCH TO APRIL CHANGES IN CIRCULATION 
AND WEATHER 
CIRCULATION 


Principal circulation changes from March to April can 
best be interpreted from figure 3, the 30-day mean 700-mb. 
anomalous height change from March to April. (The 
reader is referred to figure 4 of the preceding report of this 
series [1] for the March contour pattern.) 

Excepting Asia, the overall decrease in DN at high 
latitudes, combined with the increased DN at middle and 
low latitudes, indicates: (1) a general increase in westerly 
winds (a rising index), and (2) a weakening of the blocking 
in the Western Hemisphere. 

Blocking on the Asiatic side of the Pole in the Laptev 
Sea was transferred to the south of Novaya Zemlya in 
April, where a 270-foot positive anomalous change was 
centered. The increases over Manchuria were possibly 
a response to the upstream deepening of the mean Low 
over north-central Siberia. 

Large height changes occurred in the Pacific from 
March to April. A ridge in March, extending from the 
Gulf of Alaska to Midway Island, was replaced in April by 
a flattening of the westerlies and cyclonic curvature. An 
extensive increase in DN in the eastern Pacific resulted. 
Consequently, the west coast trough of March was forced 
inland to the west-central United States where we see a 
decrease of DN extending northward to a 580-foot maxi- 
mum change over Hudson Bay. This does not mean that 
blocking left that part of Canada completely, but it must 
have been eroded substantially. 

In the Atlantic and western Europe the increased DN 
indicates that there was a great growth of the subtropical 
ridge. This, in turn, was associated with weakening and 
retrograde motion of the western Atlantic trough to the 
west of Greenland and along the United States east coast. 
Thus it became a part of the North American trough- 
ridge complex and consequently a factor in United States 
weather. 


WEATHER 


Changes in weather over the United States from March 
to April are most easily expressed in terms of temperature 
and precipitation classes. The reports from one hundred 
selected stations used routinely are summarized in table 2 
in terms of class changes. (For example: a class change 


TaBLE 2.—Temperature and precipitation class changes from March 
to April 1958 in the United States (based on 100 selected stations). 


Temperature Precipitation 
Class change Percent of Class change Percent of 
stations Stations 
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of +2 in temperature means that a particular station 
may have warmed in the mean from normal to much above 
or from much below to normal.) Table 2 shows that 
April was considerably warmer than March. Only 7 per- 
cent of the country had lower anomalous temperatures 
during April, while 55 percent of the stations were warmer 
than in March. The remaining 38 percent showed no 
change. | 

The class changes of precipitation from March to April, 
expressed in terms of three classes in table 2, show that 
April was somewhat wetter than March. Forty-five per- 
cent of the cities showed an increase of one or two classes, 
while 21 percent totaled less precipitation than in March. 

The changes in circulation and weather described above 
did not proceed in a gradual manner during April. It is 
true that if one is given monthly mean charts at sea level 
and aloft, certain qualitative comments and deductions 
can bemade. However, there was considerable variability 
during April, and therefore, many of the weather details 
can be more easily recognized and discussed by use of 
shorter-period time-means, 


4. VARIABILITY WITHIN THE MONTH IN AND NEAR 
NORTH AMERICA 


APRIL 1-15 


The 700-mb. circulation for the first 15 days of April 
shown in figure 4a presents a picture of blocking over 
North America. The split in the westerlies approaching 
the continent not only provided for two storm tracks, but 
also left between the two belts an area of shattered trough 
and blocking ridge. 

The anomalous temperatures for the United States are 
shown in figure 5a. Here we see a pattern which is not 
atypical of a blocking regime and which compares well 
with the DN flow in figure 4a. The DN flow derived from 
the blocking was easterly across the northern States, with 
heights above normal from North Dakota through the 
northern Great Lakes. Even though some of the heights 
were slightly below normal, that type of flow encourages 
warmer than normal temperatures since in April the cold 
source region is not in eastern but in western Canada. 

The below normal temperatures in the southern two- 
thirds of the United States resulted almost. exclusively 
from intrusions of relatively cold Pacific air in a moder- 
ately fast, low-latitude, carrying current along the mean 
jet axis (fig. 2). There were only two penetrations of 
continental polar air into the country in this period, and 
those were of short duration. (See Chart IX.) Most of 
the polar anticyclonic activity was concentrated in the 
Hudson Bay area, where Highs are frequent in April [5], 
and where they became the sea level counterpart of upper- 
level blocking. 

A comparison of the 2 weeks comprising the first half 
of April? brings out some interesting occurrences. An 


2 It is not implied that the two 5-day means will total the 15-day mean shown. For 
convenience and economy of space only 4b and 4c were selected as representative of the 
first and second weeks of the month. 
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Figure 4.—(a) 15-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for April 
1-15, 1958. The blocking ridge in Canada is apparent from both 
700-mb. contours and height departures, but on the monthly 
700-mb. chart (fig. 1) only height departures showed the block- 
ing. (b) 5-day mean 700-mb. contours (tens of feet) for April 
1-5, 1958, representative of the week ending April 6. Note 
blocking High north of Lake Huron and shearing of United 
States west coast trough. (c) 5-day mean 700-mb. contours 
(tens of feet) for April 8-12, 1958, representative of week ending 
April 13. Well-defined southwestern trough occurred with ret- 
rogression of blocking and simultaneous growth of ridge into the 
Pacific Northwest. 
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TEMPERATURE ANOMALY 
a APRAL 1-15, 1958 


Midnight, Let, April 13, 1088 


Figure 5.—Departure of average surface temperatures from nor- 
mal (F. °) during the first half of April 1958. (a) April 1-15. 
This type of zonal pattern, warm in the north and cold in the 
south, represents the influence of blocking. (b) Week ending 
April 6. Cold air covered the West, but the central half of the 
country was warm as the west coast trough sheared and moved 
inland. (c) Week ending April 13. Cold air progressed to the 
east coast through the trough, and the west coast warmed as the 
700-mb. ridge extended into the Pacific Northwest. (b and ¢ 
from Weekly Weather and Crop Bulletin, National Summary, vol. 
XLV, Nos. 14 and 15, April 7 and 14, 1958.) 
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Based on preliminary telegraphic reports 


Ficure 6.—Total precipitation, inches, for: (a) April 1-6, 1958. 
Amounts over 6 inches in the San Francisco region caused local 
flooding before the trough moved inland. (b) April 8-13, 1958. 
Northerly 5-day mean 700-mb. flow over the West inhibited rain- 
fall as the trough moved inland. Cyclogenesis in the Gulf of 
Mexico early in the week was responsible for the heavy accumu. 
lation in the Southeast. (From Weekly Weather and Crop 
Bulletin, National Summary, vol. XLV, Nos. 14 and 15, April 7 
and 14, 1958.) 


inspection of figure 4b reveals a blocking High betwéen 
James Bay and Lake Superior, in phase with a tropical 
ridge in the eastern Gulf States. The northern portion 
of the west coast trough stayed well offshore, while that 
segment below 40° N. sheared and moved rapidly inland. 

Successive invasions of fairly cold maritime Pacific air- 
masses followed, and most of the West became much 
colder than normal, with temperatures averaging more 
than 9 F.° below normal over southern Nevada (fig. 5b). 
Meanwhile, the balance of the country was normal or 
above, except for below normal temperatures in the East 
which lagged the offshore trough. 

Precipitation (fig. 6a) was most notable in California. 
Record amounts (table 1) accumulating in this first week 
of the month were accompanied by severe local flooding 
and mudslides and heavy snows in the mountains [6]. 
Only light amounts of rain fell near the trough in the 
southern Plains owing to the foehn drying effect and the 
absence of Gulf moisture. 
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The second week of April was quite different from the 
first. The remnants of the offshore west coast trough 
were replaced by strong ridging (fig. 4c) into the Pacific 
Northwest as the blocking High retrograded to North 
Dakota. Thus, as the western trough-ridge couplet 
became appreciably more meridional, a temperature 
reversal took place in the second week (fig. 5c). 

As rapid cooling spread eastward from the Southwest, 
above normal temperatures in the South and Central 
States dropped to below normal. Temperatures as much 
as 14 F.° below normal were reported in the Texas Pan- 
handle, a change of some 15F.° in the mean in a week’s 
time. In the Pacific Coast States temperatures were 
above normal under the ridging aloft. 

Precipitation in the West was generally less than one- 
half an inch (fig. 6b), with an inhibiting northerly flow 
aloft to the rear of the Southwest trough. As westerlies 
with broad cyclonic curvature swept over the Gulf States, 
rainfall amounts exceeding 7 inches were deposited in 
Florida with the development of two Gulf storms. (See 
Chart X.) 

In general, then, as the week progressed, precipitation 
was most marked in the eastern and southern parts of 
the country, and average temperature anomalies assumed 
a zonal nature with warm air to the north and cold air to 
the south. Sudden changes followed, however, to give 
a sharply different picture of those elements. 


APRIL 16-30 


Changes in 700-mb. circulation during the last half of 
April were quite abrupt. Even a cursory comparison of 
the last 15 days of April (fig. 7a) with the first 15 days 
(fig. 4a) shows certain obvious gross differences in and 
near North America. 

The amplification in the eastern Pacific was perhaps 
the most important areal change in the Western Hemi- 
sphere. In the Gulf of Alaska there were rises of some 800 
feet from the previous 15-day period. That imposing 
block was probably related to the suppression of the High 
in the Arctic Ocean and the discontinuous retrogression 
of the Canadian block. Regardless of its evolution, how- 
ever, its effect downstream was of great importance. One 
result was the consolidation of the multiple trough struc- 
ture in the western United States of the first half of the 
month to an almost full-latitude trough in the central 
United States during the following 15 days. The change 
was reflected somewhat in the DN distribution (fig. 7a), 
but much more so in the changes in anomalous tem- 
peratures. 

Temperature departures from normal for the last half 
of the month (fig. 8a) hardly resemble those of the first 
half of April. Significant warming in the Southwest and 
East, coupled with vigorous cooling in the Northwest, 
produced a pattern of anomalies which was dominantly 
meridional, in contrast to the zonal pattern of April 1-15. 

The blocking ridge in the Gulf of Alaska and the accom- 
panying northerly DN flow over western North America 
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Fiaure 7.—(a) 15-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for April 
16-30, 1958. Well-defined United States trough replaced shat- 
tered troughs and blocking ridge of first half of month. Note the 
more meridional pattern of height departures from normal. (b) 
5-day mean 700-mb. contours (tens of feet) for April 15-19, 1958, 
representative of week ending April 20. Compared with the 
previous week (fig. 5c), the westerlies strengthened considerably 
in higher latitudes over North America and weakened in lower 
latitudes. (c) 5-day mean 700-mb. contours (tens of feet), for 
April 22-26, 1958, representative of week ending April 27. This 
is an example of the extreme variability within the month, with 
change from a zone! to a meridional pattern in one week. 
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Figure 8.—Departure of average surface temperatures from nor- 
mal (F. °) during the last half of April 1958. (a) April 16-30. 
This meridional distribution of temperature anomalies is in 
marked contrast to that of the previous 15-day period. (b) 
Week ending April 20. Normal or above normal temperatures 
over most of the country were a response to the northerly migra- 
tion of westerlies during this week. (c) Week ending April 27. 
Here again is another reversal. During this week cold continen- 
tal polar and maritime Pacific airmasses dominated the country 
as a result of the large amplification in the eastern Pacific Ocean. 


138 Aprit 1958 
im 
thee / for the Week Ending 
Vv 4 -Departere of Awonage from Norma! for the Week Ending 
dee 
a”, 
pla 
sou 
| air 
See] 
per 
tior 
In 
: 40° 
of 1 
sho 
cou 
tur 


05 to} Inch 


Based on preliminary telegraphic reports 


Total Precipitation, Inches, for the Week Ending 


Rave 

= 

OS ince Over 2 
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FicurE 9.—Total precipitation, in inches, for: (a) April 15-21, 
1958. Most of the precipitation was accounted for by one storm 
in the Gulf of Mexico early in the week and several frontal pas- 
sages into the Pacific Northwest. (b) April 22-28, 1958. The 
pattern of areas of precipitation exceeding 1 inch has a fair cor- 
respondence with the areas of maximum cyclonic curvature on 
the 5-day mean chart for the same week (fig. 7c). 


were basically responsible for the fairly strong advection 
of cold air into the northwestern third of the country. 
There must have been some dynamic cooling as the trough 
deepened, but this was probably of less significance than 
the advective cooling. Decreases in anomalous tempera- 
tures in the Northwest ranged as high as 16.3 F.° at 
Winnipeg, with changes in the 5°-10° category common- 
place. The eastern third of the United States was domi- 
nated by widespread warming which was caused by the 
southerly DN flow and the advection of maritime tropical 
tir in advance of the mean trough. 

Further details of the circulation of April 16-30 can be 
seen by an examination of the weekly variability in that 
period. Figure 7b shows the 5-day mean 700-mb. circula- 
tion centered in the middle of the week ending April 21. 
In this instance the westerlies were quite strong north of 
40° N., but rather flat and ill-defined in lower latitudes 
of the United States. The temperature pattern (fig. 8b) 
showed a remarkable reversal. About three-fourths of the 
country was flooded by normal and above normal tempera- 
tures with readings as high as 91° F. in South Dakota, 
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while cooler air began to invade Washington. Residual 
cold air along the Gulf of Mexico kept temperatures there 
below normal, although less so than in the previous week. 

Along with the general warming came a large area of 
light precipitation (fig. 9a). Heavy amounts were confined 
to the Pacific Northwest coast, following several frontal 
passages, and to the Southeast as a result of one Gulf storm 
early in the week. 

Figure 7c is representative of the circulation for the last 
week of April. The change from the previous week’s 
circulation was quite extreme. From the eastern Pacific 
to the central United States a strong meridional flow was 
observed, in contrast to the strong zonal flow of the previ- 
ous week, as amplification in the eastern Pacific and west- 
ern North America reached its maximum. 

Explosive deepening in the Great Basin on the 22d and 
23d was associated barotropically with the propagation of 
vorticity from the eastern Pacific and baroclinically with 
the strong thermal advection (for more details see [7]). 
It was followed by the coldest intrusion of Canadian air of 
the month. Some temperatures in the West fell as much 
as 12 F.° below normal (fig. 8c), where they had been as 
much as 9 F.° above normal the week before (fig. 8b). 
An example of the temperature gradient accompanying 
that storm, as reported by the Weekly Weather and Crop 
Bulletin [6], was the record high temperatures of 91° F. 
in eastern Kansas on the 23d, and a record minimum of 
21° F. in northwestern Kansas on the 24th. Not as 
extreme, but indicative of the weekly variability, was the 
low of 31° F. at Roswell, N. Mex., on April 14 followed by 
a maximum of 91° F. on April 22 at the same city. 

Cold air surging eastward as far as West Virginia toppled 
various minimum temperature records, while normal and 
above average temperatures were recorded in the East 
Coast States and in portions of the South and Southwest. 

The week’s precipitation is shown in figure 9b. Heavier 
accumulations were confined to the northern Rocky 
Mountain States and in a corridor extending from northern 
Texas to New England, while the Southwest was mostly 
dry. Rising flood stages were noted in parts of New Eng- 
land and the South, but serious flooding was not of major 
proportions this week. 

The variability of temperature (and to some extent of 
precipitation) during the month was intimately related 
to the variability of the 700-mb. mean flow. The relation- 
ship of another group of phenomena, severe storms, will be 
considered next. 


TORNADO ACTIVITY 


The incidence of tornadoes in April 1958 (table 3) was 
significantly lower than that reported in April 1957. 
Since acceptable normals are not available, one cannot 
classify this month’s tornado totals relative to normal. 
However, it may be of interest to briefly compare this 
year’s April circulation (fig. 1) with that of April 1957 
(fig. 8 of [8]). The positions of the central United States 
troughs were comparable, but the similarity ends there. 


Total Precipitation, Inches, for the Week Ending 
Midnight, | April 20, 1958 = 
— 
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TasLe 3.—Number of tornadoes and funnel clouds reported in the 
United States, January-April 1957 and 1958.* 


1957 1958 
Month 
Torna- | Funnels | Total Torna- | Funnels| Total 
does does 
a 18 1 19 14 4 18 
5 0 5 22 5 27 
39 26 65 17 10 27 
217 164 381 **85 **178 
Wetcededsesion 279 191 470 138 112 250 


p pee tabulated by the U. 8. Weather Bureau, Office of Climatology, Washington 
**Preliminary totals. 


In April 1957 a ridge some 200 feet above normal lay 
along the east coast. That ridge was accompanied by a 
strong southerly positive DN flow over the Midwest, a 
condition conducive to a mean influx of warm, moist, 
tropical air. The same area in April 1958 was dominated 
by below normal heights and did not have a southerly 


DN flow suggestive of extended intrusions of maritime. 


tropical air into the Midwest. 

It seems fairly clear that this year’s decrease in tornado 
activity was effected by the persistent block in eastern 
Canada, a real deterrent to middle- and low-latitude 
ridging in the eastern United States. 
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RAPID CYCLOGENESIS OVER THE GREAT BASIN, APRIL 22-23, 1958 


ROBERT F. SHAW AND ELLIS J. JOSEPH 
National Weather Analysis Center, U. S. Weather Bureau, Washington, D. C. 


1. INTRODUCTION 


The 0000 emt surface chart for April 22, 1958, gave the 
first indications of formation of a surface Low near Wend- 
over, Utah, in the Great Basin. Just 12 hours later it 
was evident that a storm of major proportions was develop- 
ing, for central pressure of the Low dropped 17 mb. and 
the circulation about the center extended over the entire 
western third of the Nation. 

The variety of bad weather produced by this storm, 
particularly in its early stages, ranged from dust storms 
inthe southern Great Plains to record-breaking snowstorms 
in the mountainous sections of Montana and Wyoming. 
Some details of the weather conditions and of record-break- 
ing occurrences are given in section 5. The rapid develop- 
ment of the Low, however, was the most spectacular fea- 
ture of this storm and this paper is primarily concerned 
with discussion of significant factors contributing to the 
cyclogenesis. 


2. BROAD-SCALE FEATURES 


Comparison of 500-mb. charts for 1200 amt on April 20 
and April 21 (figs. 1a, 1c) shows the remarkable change in 
the upper-air flow pattern over the eastern Pacific and Gulf 
of Alaska during one 24-hour period. The sequence of (1) 
essentially simple southwest flow, (2) a split in the westerlies 
with its resultant cutoff Low, and (3) repositioning of the 
major jet axis into a northwesterly current, is not uncom- 
mon in this area, but the brief period of time required for 
the complete transition is noteworthy. We shall show how 
the requisite factors for cyclogenesis in the western United 
States, though obscure at 1200 emt on April 20, became 
apparent on the 1200 emt charis of April 21. 

Long-wave trough positions prior to April 20 were: (1) 
off the east coast of Asia, (2) near the 160° W. meridian, 
and (3) in the central United States. Average wavelength 
between troughs was 70° of longitude. A new long-wave 
trough position was established east of the Kamchatka 
Peninsula by the vigorous deepening of a short-wave 
trough which moved into that area. The resulting shorter 
wavelength between the Kamchatka trough and the trough 
at 160° W. necessitated readjustment of the long-wave 
trough positions downstream, with the Kamchatka trough 
acting as an “anchor trough.” 

Constant absolute vorticity trajectories were computed 
from points along the southwest flow on the east side of 


the Kamchatka trough. The average of these trajectories, 
when superimposed on the 500-mb. chart for 1200 emr of 
the 20th (fig. 1a), suggested the “shear and reamalgama- 
tion” process for readjustment of the long-wave trough 
positions, with a new position established over the western 
United States in approximately 48 hours. This process 
has been discussed and illustrated by Namias [1]. ‘“Over- 
shooting” of the strong southwesterly winds from the 
Kamchatka trough [2] across the sharply curved contours 
of the central Pacific ridge produced an eastward shearing 
effect on the northern portion of the trough at 160° W. 
The northern branch of the westerlies consequently in- 
creased and cut off the lower portion of the trough. The 
southern branch of the westerlies weakened and eventually 
became restricted to circulation around the newly-formed 
cutoff Low. The 300-mb. jet-stream axes as shown in 
figures 1 and 2 clearly depict the changes in relative 
strength and position of the two branches of westerly flow. 

The second half of the process, “reamalgamation,”’ fol- 
lowed the breakthrough of the northern branch of the 
westerlies. The ridge over the Aleutians (fig. lb) moved 
into the Gulf of Alaska in response to the advection of 
warmer air and anticylonic vorticity into this region, while 
similar advection occurred in southeasterly current which 
developed on the northeastern side of the cutoff Low (fig. 
1d). Thus the two ridges, one from the west and the other 
from the south, combined to produce strong anticyclo- 
genesis over the Gulf of Alaska. 

With the elimination of the long-wave trough position 
at 160° W. through the process discussed above, the 
resulting wavelength between the Kamchatka long-wave 
trough and the long-wave trough in the central United 
States was in excess of the stationary value required by 
the Rossby wave formula, suggesting that a major read- 
justment of waves was necessary. This readjustment of 
the super-stationary wavelength was effected through 
what Cressman [3] terms “‘discontinuous retrogression.” 
Comparison of the 0000 emr 500-mb. charts for April 21, 
22, and 23 (figs. 1b, 1d, and 2b) shows how the long-wave 
trough position, which was over the central United States, 
retrograded to the Plateau Region. The anticyclogenesis 
in the Gulf of Alaska, combined with the necessity for 
retrogression of the long-wave trough position from the 
central United States, provided a field of motion favorable 
for the plunging of a cyclonic vorticity maximum and 
associated pool of Arctic air from Alaska southeastward 
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Figure 1.—500-mb. charts. 


(a) 1200 amr, April 20, (b) 0000 emt, April 21, (c) 1200 amr, April 21, and (d) 0000 emt, April 22, 1958. 


Contours (solid lines) at intervals of 200 geopotential feet; isotherms (dashed lines) at intervals of 5° C.; and areas of absolute vor- 


ticity (shaded) greater than 150 10-® sec.-1. 
constant absolute vorticity trajectory. 


into the Great Basin. In its broad-scale evolution, this 
situation was strikingly similar to the outbreak of Arctic 
air into the Pacific Northwest during November 1955 


[4], [5]. 
3. VORTICITY CONSIDERATIONS 


To show more clearly the source and development of the 
500-mb. trough which preceded the surface development, 
the movement of areas of maximum cyclonic absolute 
vorticity at 500 mb. were tracked (fig. 3). Absolute vor- 
ticity charts were prepared by the Joint Numerical Weath- 
er Prediction Section of the National Meteorological 
Center for 12-hour intervals from 1200 emt, April 20 


300-mb. jets denoted by heavy arrows. 


Hatched arrow in (a) represents average 


through 0000 emt, April 24. On the 500-mb. charts in 
figures 1 and 2, areas of absolute vorticity in excess of 
150 x 10~* sec. ~! are shaded. The distribution of these 
vorticity maxima in the 500-mb. flow pattern shows quite 
plainly the movements and tendencies of the upper-air 
troughs associated with the surface cyclogenesis. 

The apparent effect of the cyclonic vorticity maximum 
which moved northeastward out of the Pacific trough’ 
was temporary retardation and flattening of the ridge 
building off the west coast. The short-wave trough i 
Alaska on April 21 moved southeastward into the portion 


1 In the following discussion the trough néar 160° W. longitude will he referred to 48 
the Pacific trough. 
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Ficure 2.—500-mb. charts. (a) 1200 amr, April 22, 1958. (b) 0000 amr, April 23, (c) 1200 emt, April 23, and (d) 0000 amr, April 24. 
300-mb. jets denoted by heavy arrows and areas of maximum absolute vorticity shaded. 
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Ficure 3.—Tracks of centers of areas of maximum absolute vorticity (dashed lines), 1200 emt, April 20 to 1200 cmt, April 23, 1958, and 
track of surface Low (solid line), 0000 cmt, April 22 to 0600 amr, April 24, 1958. 
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Ficure 4.—Surface synoptic charts with 1000—-500-mb. thickness lines (dashed). 


(a) 1200 er, April 20, (b) 0000 cmv, April 21, (c) 120 


emt, April 21, and (d) 0000 Gmv, April 22, 1958. 


of the ridge which had been weakened by advection of 
cyclonic vorticity from the Pacific trough. Thus a strong 
concentration of cyclonic vorticity was propagated off the 
coast of Washington on April 22 (fig. 2a), just prior to the 
spectacular surface development downstream during the 
following 24 hours. 


4. SURFACE DEVELOPMENT 


For a week or more prior to April 20, zonal flow in the 
upper troposphere steered surface systems from the Pacific 
into the North American continent in the vicinity of the 
United States-Canadian border with some semblance of 
regularity. With the abrupt change in the upper-level 
flow from zonal to meridional on April 21, the surface 
picture also changed quite markedly (fig. 4). The pro- 


gression of surface Lows across the west coast was com- 
pletely eliminated by the block in the Gulf of Alaska, 
which was reflected at the surface by the 1040-mb. High 
in the vicinity of Ship “P” (50° N., 145° W.) on April 22 
(fig. 4d). A lobe of this High, extending southward along 
a northwest-southeast axis, dominated the surface picture 
for the remainder of the month in the area off the west 
coast of North America. Inland, the surface chart pre 
sented a rather nondescript collection of small Highs and 
Lows in the western United States. 


As the 500-mb. trough, and its associated vorticity 
maximum, moved southeastward from Alaska, surface 
pressures began falling in the Pacific Northwest about 
1800 amr on the 21st. During the following 6 hours the 
isallobaric fall center moved southeastward to the Greal 
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Ficure 5.—Sectional surface synoptic charts. (a) 0300 emt, April 22, (b) 0600 emt, April 22, and (c) 0900 amr, April 22, 1958. 
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ut Ficure 6.—Surface synoptic charts with 1000—500-mb. thickness lines (dashed), (a) 1200 amr, April 22, (b) 0000 emt, April 23, (c) 1200 
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Basin and 3-hourly pressure falls of 3 to 4 mb. were re- 
ported from most stations in this region. The first closed 
circulation of the new storm, although it probably was 
not recognized as such at the time, appeared on the 0000 
emt surface chart for April 22 just south of Wendover, 
Utah (fig. 4d). Central pressure at this time was 1009 mb. 
The series of sectional surface charts for intermediate 
synoptic hours (fig. 5) shows development during the 
next 12 hours. By 1200 amr, April 22 (fig. 6a), the Low 
had moved eastward into southwestern Wyoming and 
deepened to a 992-mb. center, a remarkable fall of 17 mb. 
in just 12 hours. 

Meanwhile, a weak cold front through central Cali- 
fornia, a remnant of a frontal system which passed through 
the Plateau during the preceding 36 hours, became revital- 
ized and extended into the developing Low (fig. 6a). The 
1,000-500-mb. thickness gradient to the north of the cold 
front increased markedly in response to the strong cold 
air advection associated with the deepening upper trough. 
The sharp discontinuity of the thickness gradient across 
the front is shown very well in figures 6a-d. 

The storm attained its greatest development in eastern 

Jolorado at 0000 amr, April 23 (fig. 6b), with a central 
pressure of 982 mb. During the 24-hour period after this 
time the Low became complex, when several small pertur- 
bations developed along the front as reflections of the dis- 
organized upper-level features (possibly complicated by 
irregular terrain) depicted in figure 2c. By 0000 emr on 
the 24th, however, both the upper-level and surface Lows 
had become organized into single vortices and had begun 
to curve northeastward toward the Great Lakes. 


5. ASSOCIATED WEATHER 


Some features of the weather associated with this storm 
are briefly described in [6]. The influx of cold air behind 
the Low as it moved northeastward brought freezing tem- 
peratures to much of the Rocky Mountain and Great 
Plains States on the 23d. A good example of the sharp 
thermal contrast across the cold front occurred in Kansas 
where maximum temperatures in the 90’s (° F.) were re- 
ported ahead of the front in the eastern part of the State 
while Goodland, in northwestern Kansas, reported a new 
record late season minimum of 21° F. 

Precipitation accompanying the storm was both wide- 
spread and varied. Rain and drizzle fell ahead of the 
warm front and thunderstorms with hail, freezing rain, 
and snow were reported behind the cold front. Strong 
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TABLE 1.— Weather occurrences with storm of 22-23, 1958, that 
equalled or broke prior records. 


Place Date Weather occurrence 
22-23 1.38 in. Greatest 24-hour tation 
amount since 1944. wus 
ee a 22-23 | Heaviest snowstorm of record so late ip 
Winslow, Ariz 23 299 recorded 
at this station. 
Salt Lake City, Utah..............- 23 15in. Heaviest 24-hour April snowfalj 
on record. 
or 
Goodland, 24 F. Record late season low tempera 
ure. 
equalled record established in ne 


winds in the warm sector raised dust storms from Arizona 
to Oklahoma. South-central Montana and the Big Horn 
Mountain region of Wyoming had one of the heaviest 
snowstorms on record April 23 and 24. Up to 3 feet 
of new snow was reported from the Big Horn Mountains, 
near Sheridan, Wyo., and in Montana, Red Lodge received 
55 inches, and Nye and Mouat Mine 60 inches. 

Table 1 lists weather occurrences with this storm that 
equalled or broke prior records. 
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Publication of Weather Notes 


Many years ago the Monthly Weather Review published detailed eyewitness 
accounts of exceptional storms. These accounts both enrich the meteorolo- 
gist’s knowledge of storms and provide him with particular details that cannot 
be found elsewhere. Because such information bears directly upon questions 
the meteorologist must attempt to answer about weather phenomena (for ex- 
ample, the identification of storms as tornadoes), and because the information 
has potential value in both the research and service programs of the Weather 
Bureau, publication of eyewitness accounts and brief analyses of exceptional 
storms and other meteorological phenomena was resumed in the April 1955 
issue. They appear from time to time under the heading ‘‘Weather Notes.” 

Contributions to these ‘‘Notes” are invited from readers of the Review. 
There is no limitation placed on length of description but it is expected that 
most will be short accounts. Any weather peculiarities, whether storms or 
other phenomena, are acceptable subject matter. Material should be 
addressed to Editor, Monthly Weather Review, U. S. Weather Bureau, 
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Washington 25, D. C. 


Description of Charts 


CHART I. A. Average Temperature (° F.) at Surface. 
B. Departure of Average Temperature from Normal.—The 
average monthly temperature presented in Chart I-A is 
computed from the average daily maximum and the 
average daily minimum which in turn are computed from 
the daily maximum and minimum temperatures reported 
by some 225 first-order Weather Bureau stations and 700 
cooperative stations. The departures from normal are 
presented in Chart I-B. They are based on the 30-year 
normals (1921-50) for the first-order Weather Bureau 
stations and on means of 25 years or more (mostly 1931-55) 
for the cooperative stations. 

CHART II. Total Precipitation.— 

CHART III. A. Departure of Precipitation from 
Normal (inches). B. Percentage of Normal Precipitation.— 
Chart II is based on daily precipitation records at about 
800 Weather Bureau and cooperative stations. In Chart 
III the anomaly in the month’s precipitation is shown as 
a departure from the normal total and as a percentage of 
the normal total. These anomalies show the deviations 
from the 30-year normals (1921-50) for about 225 first- 
order Weather Bureau stations in Charts III A and B, 
supplemented in Chart III-A by the deviation from means 
of 25 years or more (mostly 1931-55) for about 700 
cooperative stations. 


CHART IV. Total Snowfall.— 

CHART V. A. Percentage of Normal Snowfall. B. 
Depth of Snow on Ground.—Chart IV gives the total depth 
in inches of unmelted snowfall as reported during the 
month by Weather Bureau and cooperative stations. 
This is converted in Chart V-A into a percentage of the 
normal total amount computed for each Weather Bureau 
station having at least 10 years of record. The depth of 
snow on ground is that reported by both Weather Bureau 
and cooperative stations as of 7:00 a. m. Est on the last 
Monday of the month. This is reported only for the 
months December through April. The snowfall charts 
are presented each month November through April. 

CHART VI. A. Percentage of Sky Cover Between Sun- 
rise and Sunset. B. Percentage of Normal Sky Cover 
Between Sunrise and Sunset.—These charts are based on 
visual observations made hourly at Weather Bureau 
stations and averaged for the month. Sky cover includes, 
in addition to cloudiness, obscuration of the sky by fog, 
smoke, etc. Normal amount of sky cover is computed 
for stations having at least 10 years of record. 

CHART VII. A. Percentage of Possible Sunshine. B. 
Percentage of Normal Sunshine.—Chart VII-A shows the 
amount of sunshine received in terms of percentage of the 
total hours of sunshine possible during the month. In 
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Chart VII-B this is shown as a percentage of the normal 
number of hours of sunshine received; normals are com- 
puted for Weather Bureau stations having at least 10 
years of record. 

CHART VIII. Average Daily Values of Solar Radia- 
tion, Direct and Diffuse—Plotted on the chart are the 
monthly means of daily total solar radiation, both direct 
and diffuse, in langleys (gm. cal. cm.~') for all Weather 
Bureau stations which record this element. Supplemen- 
tary data, for which limits of accuracy are wider than for 
those data shown, are drawn upon in making the analysis. 
The inset shows the percentages of the mean based on the 
period 1951-55. 

CHART IX. Tracks of Centers of Anticyclones at Sea 
Level.— 

CHART X. Tracks of Centers of Cyclones at Sea 
Level.—Centers which can be identified for 24 hours or more 
are tracked in these charts. Semi-permanent features such 
as the Great Basin and Pacific Highs and Colorado and 
Mexico Lows are not shown. The 7:00 a.m. Est positions 
are shown by open circles, with the intermediate positions 
at 6-hour intervals shown by solid dots. The date is given 
above the circle and the central pressure to whole millibars 
below. A dashed track indicates a regeneration rather 
than actual movement to the next position. Solid square 
indicates position of stationary center for period shown 
beside it. 
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CHART XI. Average Sea Level Pressure (mb.) and 
Surface Windroses—The average monthly sea level 
pressure is obtained from the averages of the 7:00 a. m. 
and 7:00 p.m. est pressures reported at Weather Bureay 
stations. Windroses are based on the hourly wind 
directions (to 16 points of the compass) reported by 
Weather Bureau stations, each circle or arc indicating 
5 percent of the time. The inset shows the departure of 
the average pressure from the normal average computed 
for each station having at least 10 years record and for 
each 10° intersection in a diamond grid over the oceans 
from interpolated values read from the Historical Weather 
Maps for the 20 years of best coverage prior to 1940. 

CHARTS XII-XVII. Average Height, Temperature, 
and Resultant Winds, 850, 700, 500, 300, 200, and 100 mb.— 
Height is given in geopotential meters and temperature 
in degrees Celsius. These are the averages of the 1200 
Gat radiosonde reports. Wind speeds are given in knots; 
flag represents 50 knots, full feather 10 knots, and half 
feather 5 knots. Directions are shown to 360° of the 
compass. Winds are based on rawins at 1200 emr. 

NOTE. Tabulations of exact values of most of these 
charted elements for Weather Bureau stations are printed 
each month in Climatological Data—National Summary, 
and annual averages are presented in the Annual Issue of 
that publication each year. 
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Chart I. A. 


Average Temperature (°F.) at Surface, April 1958. 
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B. Departure of Average Temperature from Normal (°F.), April 1958. 
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A. Basedon reports from over 900 Weather Bureau and cooperative stations. The monthly average is half the sum of the monthiy 
average maximum and monthly average minimum, which are the average of the daily maxima and daily minima, respectively. 
B. Departures from normal are based on the 30-yr. normals (1921-50) for Weather Bureau stations and on means of 
25 years or more (mostly 1931-55) for cooperative stations. 
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Chart III. A. Departure of Precipitation from Normal (Inches), April 1958. 


B. Percentage of Normal Precipitation, April 1958. 
| 
| 
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Normal monthly precipitation amounts are computed from the records for 1921-50 for Weather Bureau stations and from 
records of 25 years or more (mostly 1931-55) for cooperative stations. 
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Chart V. A. Percentage of Norma! Snowfall, April 1958. 
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A. Amount of normal monthly snowfall is computed for Weather Bureau stations having at least 10 years of record. 
B. Shows depth currently on ground at 7:00 a.m. E.S.T., of the Monday nearest the end of the month. It is based on reports 
from Weather Bureau and cooperative stations. Dashed line shows greatest southern extent of snowcover during month. 


APRIL 1958 


Chart VI. A. Percentage of Sky Cover Between Sunrise and Sunset, April 1958. 
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A. In addition to cloudiness, sky cover includes obscuration of the sky by fog, smoke, snow, etc. Chart based on 
visual observations made hourly at Weather Bureau stations and averaged over the month. B. Computations 
of normal amount of sky cover are made for stations having at least 10 years of record. 
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Chart VII. A. Percentage of Possible Sunshine, April 1958. 


a ____B. Percentage of Normal Sunshine, April 1958. _ 
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A. Computed from total number of hours of observed sunshine in relation to total number of possible hours of 
sunshine during month. B. Normals are computed for stations having at least 10 years of record. 
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